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antenna  must  be  properly  matched.  This  thesis  presents  the 
design  of  a  matching  network  for  a  1-meter  monopole  antenna, 
operating  over  30-90  MHz  using  the  real  frequency  method 
(RFM) .  It  outlines  the  mathematical  steps  needed  to  determine 
the  equalizer  function,  which  ultimately  leads  to  the  circuit 
design.  The  goal  of  the  RFM,  given  the  real  frequency  data, 
is  to  optimize  the  Transducer  Power  Gain  (TPG) ,  and  minimize 
the  reflection  coefficient  or  power  lost  due  to  the  impedance 
mismatch.  A  complete  design  including  network  realization  is 
given.  However,  no  experimental  results  are  presented. 
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INTRODUCTION 


I . 

One  of  the  important  factors  in  antenna  performance  is  the 
input  impedance.  The  antenna  impedance  consists  of  both  real 
(resistance)  and  imaginary  (reactance)  parts.  The  reactive 
component  is  generally  unwanted  because  it  gives  rise  to 
stored  energy  in  the  near  field  of  the  antenna  [Ref.  8). 

The  input  impedance  is  primarily  determined  by  the 
geometry  and  electrical  size  of  the  antenna,  and  it  can  be 
significantly  different  from  the  impedance  of  the  generator. 
To  optimize  the  power  transfer  from  the  generator  to  the 
antenna,  it  is  necessary  to  insert  an  impedance  transformer 
between  the  two.  Ideally,  the  transformer  (or  matching 
network)  should  be  designed  to  eliminate  the  reactive 
component  of  the  antenna  impedance,  and  at  the  same  time 
provide  an  input  resistance  equal  to  that  of  the  generator. 
This  is  relatively  easy  to  accomplish  at  a  single  frequency, 
but  becomes  more  difficult  as  the  operating  frequency  of  the 
antenna  increases. 

Until  the  development  of  the  real  frequency  method  ( RFM) , 
a  broadband  matching  network  had  been  designed  by  an  analytic 
method  or  by  an  iterative  trial  and  error  procedure.  An 
analytic  method  requires  complex  and  rigorous  mathematics  even 
for  a  simple  network.  However,  RFM,  developed  by  Carlin  in 
1977,  made  designing  a  matching  network  simpler,  more  direct 
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and  less  complex.  It  does  not  assume  an  equalizer  topology, 
nor  does  it  require  an  analytic  dascription  of  the  load  (input 
impedance  of  the  antenna  in  our  case)  as  long  as  it  can  be 
obtained  by  some  means  [Ref.  1].  It  is  a  numerical  method 
that  only  requires  the  real  frequency  data  of  the  load  for  the 
frequency  band  of  interest  [Ref.  1). 

Although  several  approaches  have  been  published  for 
broadband  impedance  matching,  none  has  been  tailored  specially 
for  broadband  monopole  antennas.  However,  in  his  recent  work, 
Rao  [Ref.  10]  has  built  and  tested  a  matching  network  for  a 
loaded  monopole  antenna  at  HF  using  the  resistivity  profile 
developed  by  Wu  and  King  [Ref.  11].  Unfortunately,  that 
antenna  is  35  feet  long  and  is  not  an  acceptable  candidate  for 
manpack  or  vehicular  mount.  Therefore  a  1-meter  monopole 
antenna  was  chosen.  Like  a  dipole,  this  antenna  is  narrow 
band  and  it  is  "electrically  small"  in  the  very  high  frequency 
( VHF )  band.  However,  it  can  be  made  broadband  by  resistive 
loading  and  properly  designing  a  matching  network. 

This  thesis  describes  the  mathematical  basis  of  RFM  and 
uses  this  approach  to  design  a  matching  network  for  a  1 -meter 
monopole  antenna  operating  over  30-90  MHz. 
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II.  MATHEMATICAL  BASIS  FOR  RFM 


Figure  1  Simple  Two-port  Network 

It  can  be  seen  from  the  previous  equation  that  a  perfectly 
matched  network  will  have  a  gain  of  one.  However,  this  is  an 
unrealistic  design  that  is  not  achievable  in  practice.  Our 
goal  is  to  design  a  network  which  minimizes  p  and  maximizes 
the  power  delivered  to  the  load. 

For  a  given  load  impedance,  looking  at  the  Thevenin 
equivalent  circuit  from  the  loaded  port  [Ref.  1]  allows  us  to 
find  an  equalizer  impedance,  Zq{ w)  .  Expressing  the  TPG  in 

terms  of  load  impedance  and  equalizer  impedance, 

,  ,  4R,  (a>)  R~(u>) 

T(  u) )  - - £ - SL - -  (1) 

jZL\u>)  +  |“ 

where 

Rl  (<j)  =  load  resistance 

Rq{w)  =  equalizer  resistance 
ZL{ u>)  -  load  impedance 
Zq(o>)  -  equalizer  impedance. 
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B.  MATHEMATICAL  DEVELOPMENT 


The  fundamental  approach  of  RFM  is  its  use  of  real 
frequency  data  to  determine  an  equalizer  function.  In 
general,  the  impedance  of  the  network  can  be  complex: 
Zq(  <i>)=Rq(«)  +  jXq  (<•))/  where  Rq(w)  is  the  real  part  and  Xq((j)  is 
the  imaginary  part.  We  will  use  Rq(w)  to  find  Zq(u).  The 
key,  of  course,  is  finding  the  real  part.  In  a  step-by-step 
procedure,  the  next  two  sections  are  devoted  to  finding  the 
real  part  (resistance)  of  the  complex  impedance  function,  and 
the  following  section  derives  the  equalizer  impedance.  One 
thing  to  note  is  that  the  poles  of  the  equalizer  impedance 
must  be  in  the  negative  half  (left)  of  the  complex  frequency 
plane. 

1.  Linear  Combination  Approximation 

It  is  desired  to  design  a  broadband  equalizer  in  the 
frequency  range  <  u  <  <4,.  The  given  frequency  range  is 
first  partitioned  into  smaller  bands,  and  the  resistance  .s 
assumed  to  behave  linearly  within  each  sub-band.  Several  out- 
of-band  break  points  are  added  from  zero  frequency  to  the 
lower  frequency  <a  ,  and  one  frequency  break  point,  <4,  is  added 
beyond  u^.  The  choice  of  (4  depends  on  the  roll-off  desired. 

The  first  step  in  solving  for  an  equalizer  impedance 
is  to  obtain  a  linear  approximation  of  the  resistance,  Rq(<«>). 
The  values  of  Rq(u)  are  dependent  on  the  excursive 
resistances,  r*,  or  the  unknowns.  The  excursive  resistances 
are  the  ramp  values  between  each  of  the  break  points. 
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0  <  co  <  co  .  .  .  <  (o  ,  for  a  given  frequency  range  [Ref.  7]. 

The  number  of  unknowns  are  determined  by  the  break  points. 
For  example,  if  there  are  n  break  points,  there  are  n-1 
unknowns.  The  relationship  between  R  ( <o)  and  r  is  [Ref.  1] 


where 


N 

- r0 ak(w)  rk 


ak(o)) 


0, 


to^co 


(2) 


(3) 


and  r  =  DC  resistance. 

The  equalizer  resistance,  R,(co),  is  made  zero  for  co  > 
<4,  [Ref.  1],  From  equations  (2)  and  (3),  this  means  that 


n 


If  the  DC  resistance  value  r  is  available,  the  number  of 
unknowns  is  no  longer  n-1  but  n-2,  and  we  have 

Tn  =  -(ro  +  *  (5) 

We  have  only  considered  the  real  frequency  data  thus 
far.  However,  an  equalizer  impedance  function  has  both  even 
(real)  and  odd  (imaginary)  parts.  Since  the  resistance  is 
assumed  to  be  piecewise  linear  in  frequency,  the  reactance 
will  be  defined  in  the  same  manner  [Ref.  1].  As  in  the  case 
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of  the  resistance,  the  reactance  is  expressed  in  terms  of  the 
excursive  resistances 

*>*(«)  rk  (6) 


where  the  coefficients  b  are  obtained  from  [Ref.  3] 

<*>Jr 


bk(u)  = — - — - -  f  ln|i^jdy 

*(«*-«*-! )  J  y-« 


(7) 


They  can  be  written  in  a  closed  form  [Ref.  4]  as 

jbA(w)  ■=  - - - - -a)*!  (*+l)  log(x+l)  +  (x-1)  log|  (x-1)  |-21og(x)  ] 

where 


with  the  real  and  imaginary  parts  defined,  an  IMSL 
optimization  routine  ZXSSQ  can  be  employed  to  find  the 
excursive  resistances  required  to  produce  a  given  TPG.  The 
error  function  to  minimize  is  |T-T{  <«>)(,  where  T  is  the 
assumed  power  gain,  which  can  be  increased  until  the 
resistance  values  just  begin  to  become  negative. 

2 .  Rational  Approximation 

The  second  step  is  to  obtain  a  rational  function  which 
closely  approximates  the  piecewise  linear  curve  specified  by 
the  resistive  excursions  [Ref.  2].  This  is  done  so  that  a 
circuit  realization  of  the  equalizer  impedance  can  be 
determined  using  the  Gewertz  method  which  requires  a  ratio  of 
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polynomials.  For  convenience,  we  assume  that  the  DC 

resistance,  rD,  is  zero  in  the  subsequent  development. 

Previously  we  have  stated  that  Rq(a>)  must  be  non- 
negative  for  an  infinite  frequency  range  [Ref.  1].  This 
places  a  constraint  on  an  optimization  routine,  and 
constrained  optimization  is  difficult  to  handle.  This  is 
because  most  optimization  routines  are  written  for 
unconstrained  conditions  [Ref.  5]. 

Direct  use  of  the  unconstrained  optimization  will  lead 
to  positive  and  negative  values  of  resistances  which  are 
unacceptable.  To  get  around  this,  the  numerator  and 
denominator  polynomials  in  the  rational  function  approximation 

j?  («)  =  U0  »  Ai<*2  »...  =  A(u2) 

9  {1  +  fl^tt2  +.  .  .  +  Ba«2jl)  B(w2) 

are  expressed  in  terms  of  a  second  polynomial  of  the  form, 
Pn(©)  =  1  +  Xjcb  +  .  .  .  +  xnon.  The  denominator  polynomial, 

for  example,  can  be  written  as 

S(«2)=-|  ■  (9) 

Noting  that  Rq(0)=0  and  using  only  one  term  in  the  numerator 
polynomial,  the  rational  resistive  function  can  now  be  written 
as 
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Zq{3) 


P(S )  m  Mi  +  nx 
0(s)  102+112 


where  s=ja>,  m,  and  m2  are  the  even  parts  of  P(s)  and  Q(s) 
respectively,  and  n2  and  n2  the  odd  parts.  The  denominator 
polynomial,  Q(s),  is  related  to  B(a>:) 

B(«2)  =  B{-s2)  =  Q(S)Q(-S)  (13) 


where  Q(s)  has  all  of  its  roots  in  the  left  hand  plane  and 
Q(-s)  has  its  roots  in  the  right  hand  plane. 

We  now  solve  for  P(s),  whose  order  must  not  exceed 
that  of  Q(s).  Using  undetermined  coefficients  [Ref.  6],  we 
express  Zq(s)  as 


z  I.)  .  £(£L  .  c03°.c1a°-n...->c. 

q  Q(s)  s°+d13a~l+d23a~2+. .  .  +da 


(14) 


Equating  P(s)  term  by  term  to  (m!m2  -  n2n2)  and  solving  for 
coefficients  yields  P(s).  Reference  6  discusses  other 
procedures  for  the  solving  rational  function  of  a  driving 
point  impedance. 

4.  Circuit  Realization 

Now  that  Zq(s)  is  known,  a  circuit  that  provides  the 

required  impedance  is  obtained  by  a  conventional  synthesis 
method.  This  is  a  procedure  by  which  a  network  is  generated 
from  a  given  input/output  relationship  [Ref.  8].  The  details 
will  be  discussed  in  the  next  chapter. 
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III.  APPLICATION  OF  RFM 


In  this  chapter,  we  will  apply  the  mathematical  procedures 
of  the  previous  chapter  to  design  a  realizable  circuit.  As  an 
illustration  of  the  method,  the  results  presented  in  [Ref.  1] 
will  be  duplicated  and  then  applied  to  a  1 -meter  monopole 
antenna.  In  order  for  this  antenna  to  operate  in  a  broad 
frequency  range,  the  matching  network  must  make  the  antenna 
impedance  less  sensitive  to  frequency.  This  is  discussed 
briefly  in  the  1-meter  monopole  design  section. 

A.  EXAMPLE  OF  RFM  APPLICATION 

In  order  to  verify  a  computer  program  (Appendix  A)  and  to 
evaluate  an  IMSL  optimization  routine  (Appendix  B) ,  published 
data  generated  by  Carlin  [Ref.  1]  were  used  to  design  an 
equalizer  network. 

The  notching  network  we  wish  to  design  for  the  given  load 
is  shown  in  Figure  2.  The  normalized  frequency  range  of 
interest  is  from  0  to  1.25  (0  <  w  <  1.25),  and  an  increment  of 
0.25  will  be  used.  This  gives  6  break  points  (observation 
points).  Since  the  design  requirement  states  that  the  TPG 
must  be  maintained  at  T(0) =0.846,  this  forces  the  circuit  to 
have  a  resistance  value  of  r0=2.29  ohms.  Calculation  for  the 
DC  resistance  value  can  be  obtained  with  the  following 
formulas  [Ref.  1]: 
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Equator  Load  impedance  j 

Figure  2  Given  Circuit  with  Load 

ro  *  *j(0)  l*o  *  (15) 

°  T(oT 

where  T{0)  is  the  DC  gain. 

Based  on  these  formulas,  there  are  two  cases  of  finding 
the  rational  function:  case  1  where  rc  >  Rx ( 0 )  and  case  2 
where  r0  <  Rx(0).  In  this  example  we  will  perform  the  design 
with  case  1.  With  rc  known,  the  unknowns  are  no  longer  5  but 
4. 

Keeping  in  mind  the  concept  of  RFM,  the  load  impedance 
values  are  first  obtained  from  the  given  RLC  values  shown  in 
Figure  2.  The  results  are  given  in  Table  I.  The  ak  and  bk  of 
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equations  (3)  and  (8)  were  computer  programmed,  a  listing  of 
which  is  provided  in  Appendix  A. 


Table  I:  Impedance  of  the  Load 


Freq 

Impedance  Value  ( Q) 

0.00 

1.0000 

♦jO.OOOO 

0.25 

0.9174 

+j0.2998 

0.50 

0.7353 

+ jO . 7088 

0.75 

0.5525 

+ j 1 . 2278 

1.00 

0.4098 

+j 1 . 8082 

1.25 

0.2358 

+j3 . 0255 

1.  Linear  Combination  Approximation  of  Equalizer 
Resistance 

Substituting  equations  (2)  and  (6)  into  equation  (1), 
the  TPG  is  redefined  in  terms  of  rk  [Ref.  2] 


T(O) 


4^(0)  tr0  + 
- 


£ 


«*(»>*** 


+  zQ  *  g 


+  lx j  («)  + 


and  the  function  to  minimize  is  |r0-T{«)|.  Programming  the 

above  equation  using  an  IMSL  optimization  subroutine  ZXSSQ 
with  j  as  a  minimization  function,  the  rk  values  were 

obtained.  The  rk  values  change  with  the  initial  conditions 
provided  to  the  ZXSSQ  subroutine.  For  this  example,  the 
initial  values  were  all  set  to  zero.  These  values  in  turn 
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were  used  to  calculate  the  resistance  and  reactance  at  each 
breakpoint . 

2 .  Rational  Approximation 

Now  that  we  have  represented  Rq(<i>)  as  a  linear 
combination/  the  resistance  values  are  used  to  calculate 
Rg{u>)  .  The  function  to  minimize  is  \£g  -  at  the  discrete 

frequencies  k»0,l,.  .  .  n.  Again,  the  IMSL  subroutine 
ZXSSQ  was  used.  The  rational  function  is  obtained  as 

R  (w)  „  _ 2^29 _  =  A(tt2) 

q  1  +  4.8w2  -  10.2«4  +  8.39«®  B(w2)  ’ 

A  plot  of  the  rational  function  and  linear  combination  is 
given  in  Figure  3.  As  can  be  seen  from  the  graph,  the  piece- 
wise  linear  approximation  &  id  rational  approximation  are  in 
agreement . 

3.  Application  of  Gewertz  Method 

We  now  have  the  real  part  of  the  equalizer  impedance. 
From  the  relationship  between  £q  and  £>{s)  as  defined  in  the 

equations  (12)  through  (14),  we  can  express  B(«2)  in  terms 
of  B(-s2)  as 

B(-s2)  =  1  -  4.8s2  -  10.2s4  -  8 . 4S* . 

Finding  the  roots  of  B(-s2)  and  writing  Q(a)  in  factored  form, 
we  obtain  0(8)  as 
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C(s)  =  (s+. 887+j. 316) (S+. 887-j. 316) (s+J. 389) (s-j. 389) 
=  2.9s3+2.9s2+3.3s+l. 


The  roots  of  B{-sl)  were  found  by  a  root  finding  IMSL 
subroutine  called  PLROC. 

The  next  step  is  finding  the  coefficients  of  P(s). 
The  equalizer  impedance,  Zq{s)  ,  defined  in  terms  of  P(s)  and 

Q{s)  is 

z  (S)  =  £i£l  =  c0s3^c1s2^c2s^c3 

q  0(S)  2 . 86s3+2 . 86s2+3 . 27s+l 

Equating  the  real  part  of  Zq(s )  \0mja  to  ^{w)  ,  we  have 
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n^n^-n^  -  (c1s2+c3)  {2 . 86s2+l)  -  (c0s3+c2s)  (2 . 86s3+3 .29s)  |#.^w 
=2.29. 

Equating  the  coefficients  of  like  powers  on  both  sides,  we  get 

-2 . 86  C0S6 \a.]u  =  0  3.27ce>-2.86CjS4+2.86c1S4|,.iM*0 

CjlS2+2 . 86c3-3 . 27c2|5.^w*0  c3*2 . 29  . 

Solving  the  above  equations,  we  obtain  Zq{s)  as 

Z  (s)  =  2 . 89sa+2 ■ 89s+2 . 29 

q  "  2 . 86s3+2 . 86s2+3 . 27s+l  ' 

4.  Circuit  Realization 

From  Zq(s)  ,  it  is  necessary  to  find  the  circuit 

elements  required  to  realize  the  matching  network.  For  this 
example,  the  degree  of  the  polynomial  in  the  denominator  is 
larger  than  that  of  the  numerator.  In  order  to  divide  a 
smaller  degree  into  a  larger  degree,  we  will  convert 
Zq{s)  to  Yq(s)  , 

y  (s)  =  1  =  2  ■  86s3+2 . 86 S2 +3 . 27s-*-! 

9  ~  Zq{3)  ~  2 . 89s 2 +2 . 89s +2 . 29 

The  division  process  is  as  follows: 
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0.99s 

2 . 89s2+2 . 89S+2 . 29)2.86s3+2.86sJi+3.27s+i 

2. 86  s3 +2.86  s2 +2. 27s 

1 . 00S+1 . 


The  first  circuit  element  is  a  capacitor.  Again  converting 
and  repeating  the  process  gives 


2.89s 

1 . 0s+l)2 . 89sz+2 . 89S+2 . 29 
2 . 89s2  +2 . 895+0 . 00 
2 . 29 


and  the  second  element  is  an  inductor.  This  process  is 
continued  until  it  is  complete  and  further  division  cannot  be 
carried  out.  The  last  division  gives  us 


0.43s 

2.29)1  .OOS+I 
1 . 00s+0 
1 

which  is  another  capacitor.  The  remaining  value  is  the  DC 
resistance  which  equals  to  the  original  value  we  have 
calculated  based  on  the  assumed  TPG  of  0.846. 

Now  that  we  have  our  circuit  elements,  the  question  is 
how  these  elements  are  positioned.  The  crucial  step  is 
placing  the  first  element,  for  other  elements  follow  in  an 
alternating  sequence  of  parallel  or  series  arms  away  from  the 
load.  This  provides  a  ladder  network.  The  final  circuit  to 
achieve  Zq{s)  is  shown  in  Figure  4. 
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Load 


Equalizer 


Figure  4  Final  Matching  Network 


B.  MATCHING  OF  1-METER  MONOPOLE  ANTENNA 

We  have  gone  through  an  example  of  how  a  matching  network 
is  designed.  We  will  apply  this  procedure  to  a  l-meter 
monopole  antenna  operating  over  30-90  MHz  with  the  break 
points  chosen  at  10  MHz  increments.  The  break  points  are 
normalized  to  90  MHz.  Details  of  the  calculation  as  shown  in 
the  previous  section  will  be  avoided,  and  only  the  highlights 
will  be  presented. 

1.  Wide-banding  1-meter  Monopole  Antenna 

An  antenna  is  defined  as  broadband  "when  its  impedance 
and  pattern  do  not  change  significantly  over  about  an  octave 
or  more",  or  when  the  ratio  between  the  upper  frequency  and 
the  lower  frequency  is  greater  than  2  [Ref.  9], 

The  input  impedance  is  highly  dependent  on  the 
frequency  of  operation. 
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For  an  electrically  small 
antenna  to  operate  over  a 
broad  frequency  range  without 
continuous  fluctuation  in 
impedance  (which  in  turn 
restricts  the  power  transfer 
from  the  generator  to  the 
antenna) ,  the  antenna  must  be 
made  lossy;  i.e.,  a  resistive 
load  (or  loads)  mist  be  added 
to  the  antenna  [Ref.  9].  The 
antenna  we  wish  to  look  at  has  a  height  of  1  meter  and  a 
radius  of  0.005  meter.  It  must  operate  in  a  frequency  range 
of  30-90  MHz.  Since  the  wavelength  near  the  low  frequency  end 
is  much  larger  than  its  length,  the  antenna  is  considered 
electrically  small. 

In  order  to  calculate  the  resistive  values  to  make  the  l- 
meter  monopole  broadband,  we  used  the  concept  of  resistive 
loading  proposed  by  Wu  and  King  in  [Ref.  11].  They  used  a 
continuously  distributed  load  of  the  form 


Figure  5  Monopole  Antenna 
Divided  into  N  segment 


zHz) 


60* 

iTHTT 


(16) 


where  h  is  the  height  of  an  antenna,  and  z  is  an  incremental 
distance  from  one  end  point  of  an  antenna  (2=0)  to  the 
opposite  end  (z=h)  as  shown  in  Figure  5. 
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The  quantity  i|»  is 


ijj  =  2{sinh_1(  — )  -C(2A,2kh) -jS(2A,2kh)  +-2-  (i-e-jl2m)  (17) 
a  kh  ' 


where  a  is  the  radius,  A=ka,  and  k  is  the  wavenumber  in  free 
space  ( k  =  co,/e0n0)  .  The  quantities  C(a,x)  and  S{a,x)  of 

equation  (17)  are  defined  as  [Ref.  11) 

C(a,x)  =  f  1"c°sWdu  (18) 

o 

S (a ,x)  -  j-2i2!?du  (19) 

o 

where 

V  =  (u2  *  a2) 1/2 .  (20) 


We  have  calculated  the  various  parameters  for  a 
continuously  distributed  load  at  the  geometric  mean  frequency 
of  52  MHz.  For  the  1-meter  monopole,  ¥  is  [Ref.  13] 


«  2 


sinh"1  (200)  -  <7(0.0109,2.176)  -  ^5(0.0109,2.176) 

(l_ej<-2.176>) 


1.088 
=  9.24  -  jl. 92. 


Referring  back  to  equation  (16),  our  continuous  load  value 
using  a  30%  multiplication  factor  is  [Ref.  13] 

ZHz)  =  15  111.4  -j2.  s )  _ 
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Since  we  are  interested  in  lumped  loading,  we  obtained 
a  discrete  approximation  to  the  continuous  profile  over  a 

segment  A z,  where  N  is  the  total  number  of  segments  (  =  -k/Az) 

and  nAz  is  the  location  of  Zn.  We  have  used  N=8  for  the 
1 -meter  monopole. 

Once  the  required  load  values  were  calculated  for  each 
segment,  the  WIRE  program  [Ref.  14]  was  used  to  generate  the 
impedance  characteristics  of  the  antenna.  Table  II  shows  the 
impedance  characteristic  of  the  antenna  with  and  without  the 
load  added.  The  antenna  impedance  characteristics  plotted  on 
a  Smith  chart  are  shown  in  Figure  6 . 


Table  II:  Unloaded  and  Loaded  Impedance  for  1-meter 
Monopole  Antenna _ _ _ 


Freq  Loaded  Unloaded 

(MHz)  (ohms)  (ohms) 


(30) 

(82.57 

- j375 . 4 ) 

(3.850 

- j346 . 2 ) 

(40) 

(90.22 

- j259 .5) 

(7.040 

-j223.4) 

(50) 

(100.9 

- jl85 .0) 

(12.75 

- j!38.4) 

(60) 

(115.1 

-J132.5) 

(20.85 

-j70. 40) 

(70) 

(132.9 

—  j  94 . 32) 

(33.35 

- j9 .050) 

(80) 

(154.1 

-j67.81) 

(53.30 

+ j51 . 00 ) 

(90) 

(177.5 

-j52.43) 

(86.90 

+jll5.0) 

2.  Equalizer  Impedance  Calculation 

Again  the  same  procedure  as  above  was  used  to 
calculate  the  impedance.  Here  we  varied  the  DC  transducer 
gain  until  the  resistance  values  just  approach  zero  from 
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Figure  6  Impedance  Characteristics  of  1-meter  Monopole 

positive  values.  The  corresponding  gain  is  then  treated  as 
optimum.  The  same  program  used  to  find  the  resistance  values 
for  the  previous  example  was  used  in  this  case  after  slight 
modification.  The  difference  is  that  the  design  called  for  a 
fixed  TPG  of  0.846  for  the  previous  case,  whereas,  here  we  are 
interested  in  the  optimum  TPG.  An  optimum  TPG,  T0,  was 
located  at  0.4785,  and  the  comparison  of  the  resistance  plots 
is  given  in  Figure  7.  The  plot  shows  that  the  rational  and 
linear  approximations  closely  follow  each  other. 

The  rational  resistance  function  for  the  monopole  is 

R  (0J)  =  _ 3.590)2 _  =  A(oa) 

q  1-8 . 59 to* +56 . 64(i>4-95 . 01w6+77 . 03o»*  £(«*)  ‘ 
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Figure  7  Comparison  of  Resistance  Characteristic  of  l  meter 
Monopole  Antenna 


Following  the  steps  specified  in  the  previous  chapter,  the 

positive  real  roots  of  B(-s2)  are 

<S  +  0.3147  +  jO.  7979) 

(s  +  0.3147  -  jO.  7979) 

(S  +  0.1947  +  jO.  3420) 

(s  +  0.1947  -  J0.3420) 


which  gives  Q(s)  as 


Q(S)  =  S*  +  1.02s3  +  1.14S2  ♦  0.384S  +  0.114. 


To  be  consistent  with  the  assumed  form  of  B( o>2)  ,  £>(s)  must  be 
divided  by  a  constant  value  such  that  the  term  independent  ofs 
is  equal  to  1.  For  the  above  equation,  we  divided  by  0.114  to 
obtain  the  new  Q(s)  as 
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0(s)  =  8 .78s4+8 . 94s3+9 . 97s2+3 . 37s+l . 


The  impedance  Zq{s) 


is  now 


Assumed 


to  be  of  the  form 


Cc,S4+C1S3+C2S2+C3S+C4  _  p(s) 

~  8.78s4+8.94s3+9.97s2+3.37 s+1  ~  G(s) 


The  unknowns  in  the  numerator  can  be  obtained  by  multiplying 
the  odd  and  even  parts  and  subtracting 

J(8.78s4+9.97s2+1)  (c0s4+c2s2+c4)] 
jUj/Bj  i3i/1i*|_(8  94S3+3 ,37s)  (c1s3+c3s)  J‘ 

The  result  is  equated  to  A(w2)  for  s«ju.  Expressing  the 
unknown  coefficients  into  a  matrix  form  {AX}={B}  and  solving 
for  {X},  where  {X>  represents  the  coefficients  of  .(s),  we 
have 


1.00  0.00  0.00  0.00  0.00 

‘<?4 

0.00' 

0.00  0.00  -8.78  8.94  0.00 

c3 

0.00 

0.00  -8.94  9.97  -3.37  0.00 

<?2 

* 

0.00 

0.00  3.34  -1.00  0.00  0.00 

3.59 

0.00  0.00  0.00  0.00  8.78, 

C0, 

p.oo. 

The  IMSL  subroutine  LEQIF  is  used  to  solve  the  matrix 
equation,  the  coefficients  are  found  to  be 

Co=0.0  Ct=2.33  C2  =2 . 3  8  C3=1.77  c4=0.0 


and  zq(3) 
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_ 2.33  S3  +  2,38s2  +  1.77 _ 

8.78S4  +  8.94S3  +  9.97s2  +  3.37s  +  1 


Zg(s)  = 


In  the  above  equation,  frequency  has  been  normalized  such  that 
s=l  corresponds  to  90  MHz.  Furthermore,  the  impedance  itself 
is  normalized  such  that  Zq  =  1.0  corresponds  to  500  ohms. 

The  matching  circuit  is  now  obtained  by  the  synthesis 
method.  Since  the  power  of  the  denominator  is  larger  than  the 
numerator,  we  convert  Zg{s)  to  Yq(s)  and  reduce  the  equation 

to 


Yq{s) 


■  3.77s  + 


0.71S  + 


3.26s  + 


1 . 10s 


0.299 


To  obtain  the  value  of  the  circuit  elements,  they  need  to  be 
denormalized.  If  we  were  to  match  this  to  a  simple  75  ohm 
coaxial  transmission  line,  a  transformer  could  be  used. 
Taking  the  denormalized  resistance  value  to  be  150  ohms 
(0.299x500)  from  the  above  equation,  and  using  the  transformer 
with  a  turns  ratio  of  1:  1.41,  the  circuit  in  Figure  8  is 
obtained. 

The  TPG  with  the  matching  network  is  plotted  in  Figure 
9.  The  equalizer  tunes  the  resistance  and  the  reactance  values 
of  the  load  and  maintains  an  approximate  gain  of  0.4785. 

The  resistance  function  indirectly  determines  the 
number  of  circuit  elements  required  to  design  a  matching 
network.  From  equation  (10),  it  is  seen  that  the  maximum 


25 


Tronducer  Power  Coin 


Figure  8  Matching  Network  for  1-meter  Monopole  Antenna 
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Figure  9  Transducer  Power  Gain  ~~ 
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power  of  the  denominator  is  2  n.  The  number  of  circuit 
elements  required  for  the  design  is  usually  greater  than  or 
equal  to  n.  Therefore,  increasing  the  degree  of  denominator 
will  increase  the  number  of  circuit  elements. 

Regardless  of  what  the  highest  degree  of  the 
denominator  polynomial  is,  the  resistance  of  the  rational 
function  approximation  must  closely  follow  the  linear 
approximation.  Some  of  the  higher  order  approximations  with 
n=5  and  n=6  for  the  1-meter  monopole  are  shown  in  Figure  10. 
It  can  be  seen  from  the  figure  that  as  more  terms  are  included 
in  the  polynomial  a  closer  approximation  is  achieved; 
however,  the  TPG  is  not  significantly  affected  by  the  highest 
power  of  the  rational  resistance  function  beyond  a  certain 
number  of  terms.  Therefore,  the  minimum  acceptable  order  in 
the  rational  resistance  polynomial  should  be  used.  Otherwise, 
the  mathematics  becomes  cumbersome. 

As  an  example,  if  we  had  designed  a  matching  network 
using  n=6,  the  rational  resistance  function  would  have  ’  een 

1?  (W)* _ LZL«2 _ 

q  1-19 .  3(02 +180 . 7  to4 -73 3 . 5  to6 +1500. 8  to8 -14  2  5 .3(O10  +  505.5(O12 

and  the  equalizer  impedance  would  be 

Z  (s)  40.79ss+33.39s4+63.77s3+24.51s2  +  17,07s 

q  4.74s6+3.88s5  +  10.38s4  +  5.28s3+6.32s2+1.54s+l  * 
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Figure  10  Resistance  Curves  for  n=5  (top)  and  n=6  (bottom) 


The  above  impedance  function  is  realized  by  the  network  shown 
in  Figure  11.  The  complexity  of  the  matching  network  has  been 
increased  (compared  Figures  8  and  11) ,  but  the  TPG  has  not 
changed  significantly.  It  is  still  given  by  the  piecewise 
curve  of  Figure  9. 


28 


Figure  11  Matching  Network  for  Rational  Resistance  Function 
of  the  Highest  Power  of  12  (n*6) 
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IV.  CONCLUSION 


The  real  frequency  method  provides  an  elegant,  yet  simple 
way  of  designing  a  matching  network  for  any  load.  This  thesis 
concentrated  on  using  this  method  to  design  a  matching  network 
for  a  1-meter  monopole  antenna.  It  can  be  seen  from  the 
verification  of  Carlin's  data  and  from  the  application  to  a  l- 
meter  monopole  antenna  that  this  numerical  technique  is 
readily  realizable  and  easy  to  implement.  The  basic 
assumption  is  that  the  antenna  impedance  characteristics  are 
known  and  that  these  characteristics  are  provided  to  the 
circuit  designer.  With  the  1-meter  monopole  antenna  we  had  to 
add  a  resistive  load  along  the  antenna  in  order  to  make  its 
impedance  characteristics  broadband.  Then  we  were  able  to 
design  the  equalizer.  The  key  to  finding  the  equalizer 
impedance  is  determining  the  resistance  function,  which  is 
done  using  the  RFM.  Once  this  is  found,  the  equalizer  can  be 
completely  determined. 

The  number  of  elements  required  to  design  a  matched 
circuit  is  normally  determined  by  the  rational  resistance 
function  as  shown  in  equation  (10).  Generally,  n  indicates 
the  number  of  circuit  elements  required  for  the  design. 
Carlin  used  2n=6  for  the  highest  degree  of  the  rational 
resistance  function,  and  his  circuit  was  composed  of  three 
elements.  For  the  1-meter  monopole,  2n=8  gave  a  minimum  of 
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five  (transformer  inclusive)  circuit  elements.  For  the  higher 
orders  (n  >  6),  the  1-meter  monopole  required  seven  circuit 
elements.  The  number  of  circuit  elements  required  is  usually 
greater  than  or  equal  to  n.  Regardless  of  the  power  of  the 
rational  resistive  function,  the  TPG  is  still  the  same  as 
specified  by  the  straight  line  approximation.  Therefore,  a 
minimum  order  of  power  of  the  rational  resistance  function 
that  closely  follows  the  straight  line  approximation  should  be 
used.  Otherwise,  the  mathematics  become  cumbersome. 

For  this  thesis,  simple  software  was  written  to  calculate 
some  of  the  values.  The  ease  of  the  programming  was  due  to 
the  availability  of  the  required  software  in  the  IMSL  Math 
Library.  Listings  are  included  in  Appendix  B. 

Although  this  thesis  was  limited  to  a  simple  1-meter 
monopole  antenna,  the  techniques  presented  herein  can  be 
adapted  to  design  a  wideband  matching  network  for  any  load. 
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APPENDIX  A 

PROGRAM  FOR  1 -METER  MONOPOLE  ANTENNA 


C  THIS  PROGRAM  TAKES  THE  IMPEDANCE  DATA  OF  1  METER  MONOPOLE 

C  ANTENNA  AND  CALCULATES  THE  RESISTANCES  USING  TWO  METHODS: 

C  LINEAR  COMBINATION  AND  RATIONAL  FUNCTION.  THE  RATIONAL 

C  FUNCTION 

EXTERNAL  CURSM1 

PARAMETER  (M=11,N=9 , XJ= (N+l) *N/2 , WO=5*N+2*M+XJ) 

INTEGER  IXJAC , NSIG , MAXFN , IOPT , INFER , IER 

REAL  PARM(4) ,X(N+1) ,F(M) , XJAC (M, N+l) , XJTJ (XJ) ,WORK(WO) , 
+EPS, DELTA, SSQ,AK(20, 20) ,BK(20,20) ,Wr20) ,Z(20,20) ,Y(20,20) ,PI, 
+ AZ  (20,20) ,BZ(20,20) , CZ ( 20,20) , AY (20,20) ,  BY  (20,20)  , CY (20,20) , 
+RQ (20) , RR , XQ ( 2  0 ) , XX , T ( 2  0 ) ,RO.TPG 
COMPLEX  IMPED (20) 

COMMON  RO, IMPED, AK,BK,T,TPG 

OPEN (UNIT=1 , FILE= • RECURSM  DAT ' , STATUS= ' OLD 1 ) 

OPEN (UNIT=2 , FILE= ' FREQM  DAT  * , STATUS= ' OLD ' ) 

OPEN (UNIT=3 , FILE= ' A KM  DAT  * , STATUS= • OLD ' ) 

OPEN (UNIT=4 , FILE= ’ LOAD  DAT ' , STATUS= ' OLD ' ) 

OPEN (UNIT=8 , FILE= *  BKM  DAT ' , STATUS* ' OLD ' ) 

OPEN (UNIT=5 , FILE* ' RQM  DAT ' , STATUS* • OLD 1 ) 

READ  (6,*)  TPG 

M=ll 

N=9 

IXJAC=M 
NSIG=5 
EPS=0 . 0 
DELTA* 0 . 0 
MAXFN=2000 
IOPT=l 

DO  110  1*1, N 
110  X(I)*0.0 


R0=0 . 0 
PI=3 . 4 12 

C  ***********  CALCULATE  AK(W)  *************** 

READ (4,*) (IMPED (I) ,I=1,M) 

READ ( 2 , * ) (W(I) , 1=1, M) 

DO  10  K=1,M 
KK=K-1 

DO  20  1=1, M 

IF  (W(K)  .LE.  W(I))  THEN 
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AK(K, I) — 1. 0 

ELSEIF  (W(KK).LE.  W ( I ) . AND. W ( I )  .LE.  W(K))  THEN 
AK (K,I)=(W(I) -W ( KK ) )/(W(K)-W(KK) ) 

ELSEIF  (W(I)  .LE.  W(KK))  THEN 
AK(K, I) =0 . 0 

ELSE 

AK(K, I)=0. 0 
ENDIF 

WRITE ( 3 , * ) AK ( K ,  I ) 

20  CONTINUE 

10  CONTINUE 

C  ***************  CALCULATE  BK(W)  ****************** 

DO  50  K=1,M 
DO  51  1=1, M 
Z (K, I) —0. 0 
Y (K, I) =0 . 0 
AZ (K, I) =0 . 0 
BZ (K, I ) =0 . 0 
CZ (K, I) =0. 0 
AY (K, I) =0. 0 
BY (K,  I) =0 . 0 
CY (K, I ) =0 . 0 
BK (K, I ) =0 . 0 
51  CONTINUE 

50  CONTINUE 

DO  30  K=1 , M 
DO  31  1=1, M 

IF  (W(K)  .LT.  .001  .OR.  W(I)  .LT.  .001)  THEN 
Z (K, I) =0 . 0 
GO  TO  32 
ENDIF 

IF  (W(I)  .EQ.  W(K) )  THEN 
BZ (K, I) =0 . 0 

AZ(K,I)=(W(I)/W(K)  +1) *LOG(W(I)/W(K)  +1) 
CZ(K,I)=W(I)/W(K) *LOG(W(I)/W(K) ) 

Z (K, I) =W(K) * (AZ (K, I)  +  BZ (K, I)  -  2*CZ(K,I)) 

GO  TO  32 
ENDIF 

AZ(K,I)=(W(I)/W(K)  +  l)*LOG(W(I)/W(K)  +1) 

BZ (K, I ) = (W (I) /W(K) -1) *LOG { ABS (W ( I ) /W (K) -1) ) 
CZ(K,I)=W(I)/W(K) *LOG(W(I)/W(K) ) 

Z (K, I) =W(K) * (AZ (K, I)  +  BZ (K, I)  -  2*CZ(K,I)) 

32  KK=K-1 

IF  (KK  .EQ.  0.0)  THEN 
BK(K, I) =Z (K, I) 

WRITE ( 8 , * ) BK ( K , I ) 

GO  TO  31 
ENDIF 
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IF  (W(KK)  .LT.  .001  .OR.  W(I)  .LT.  .001)  THEN 
Y(KK,I)=0.0 

BK(K,I)=l/< (W(K)-W(KK) )*PI) * (Z (K, I) -Y (KK, I) , 
WRITE (8,*)BK(K,I) 

GO  TO  31 
ENDIF 

IF  (W(I)  .£Q.  W ( KK ) )  THEN 
BY (KK, I) =0 . 0 

AY (KK, I) = (W ( I ) /W (KK) +1 ) *LOG (W ( I ) /W ( KK) +1 ) 

CY (KK, I ) =W (I ) /W(KK) *LOG(W{I)/W(KK) ) 

Y (KK, I) =W(KK) * (AY (KK, I) +BY (KK, I ) -2*CY (KK, I ) ) 
BK(K,I)=1/((W(K)-W(KK) )  *PI) * (Z (K, I) -Y (KK,  I) ) 
WRITE ( 8  ,  * ) BK ( K ,  I ) 

GO  TO  31 
ENDIF 


AY(KK,I)=(W(I)/W(KK)+l)*LOG(W(I)/W(KK)+l) 

BY (KK, I)  = (W (I) /W(KK) -1) *LOG ( ABS (W (I ) /W (KK) -1 ) ) 
CY(KK,I)=W(I)/W(KK)*LOG(W(I)/W(KK) ) 

Y  (KK,  I)  — W (KK)  *  (AY  (KK,  I )  +BY  (KK,  I )  “2 *CY  (KK,  I )  ) 

BK(K,I)=l/( (W(K)-WTKK) ) *PI)  * (Z (K, I) -Y (KK, I) ) 

WRITE ( 8 , * ) BK ( K , I ) 

31  CONTINUE 

30  CONTINUE 

C  ************CALL  IMSL  ZXSSQ  *************** 

CALL  ZXSSQ (CURSM1,M,N ,NSIG, EPS , DELTA, MAXFN, IOPT, PARM, X, SSQ, F, 
+XJAC , IXJAC , XJTJ , WORK , INFER , IER) 

X ( 10) =- (X ( 1) +X (2 ) +X ( 3 ) +X (4 ) +X ( 5) +X ( 6 ) +X (7 ) +X ( 8 ) +X ( 9 ) ) 

WRITE (1,*) 'X(IO) ' ,X(10) 

WRITE ( 1 , *) 'T' , T 

C  **************  CALCULATE  RQ  AND  XQ  ************** 

DO  300  1=1, M 
RQ (I) =0 . 0 
300  CONTINUE 

DO  80  1=1, M 
RR=0 . 0 

DO  81  K=1 , N+l 
KK=K+1 

RR=  AK(KK, I) *X (K) +RR 

81  CONTINUE 

RQ(I)=RO+RR 
WRITE ( 5 , * ) RQ ( I ) 

WRITE (1,*) 'RQ(I) ' ,RQ(I) 

80  CONTINUE 

DO  301  1=1, M 

34 


o  o  o 


XQ ( I ) =0 . 0 
301  CONTINUE 

DO  90  1=1, M 
XX=0 . 0 

DO  91  K=1 , N+l 
KK=K+1 

XX=BK (KK, I ) *X(K)+XX 
91  CONTINUE 
XQ(I)=XX 

WRITE ( 1 , * ) 'XQ(I) ' , XQ ( I ) 

90  CONTINUE 

WRITE (1,*) 'X' ,X 
WRITE ( 1 , * ) 1 SSQ 1 , SSQ 
END 

C*** *********************  SUBROUTINE  CURSM1***** ************ ****** 


THIS  IS  A  CALLING  SUBROUTINE  TO  IMSL  SUBROUTINE  ZXSSQ.  IT 
WILL  PROVIDE  THE  LINEAR  RESISTANCE  VALUES 

SUBROUTINE  CURSM1 (R,M,N, F) 

INTEGER  M,N,I,K,N6,J 

REAL  R(N) ,F(M) ,SUM(20) ,AK(20,20) ,BK(20,20) ,T(20) ,W(20) , 

+TT (20) , RX (20) , TPG 

COMPLEX  IMPED (20) 

COMMON  RO, IMPED, AK,BK,T, TPG 

N6=N+1 
RX(N6)=RO 
DO  10  J=1,N 
RX(J)=R(J) 

10  RX (N6) =  (RX(N6)+R(J) ) 

RX(N6)=-RX(N6) 

DO  100  1=1, M 
SUMA=0 . 0 
SUMB=0. 0 

DO  15  K=1,N6 
KK=K+1 

SUM ( I) =AK (KK, I ) *RX (K) 

SUMA=SUMA+SUM ( I ) 

SUM ( I ) =BK ( KK , I ) *RX (K) 

SUMB=SUMB+SUM(I) 

15  CONTINUE 

TT ( I ) =4  *REAL ( IMPED ( I ) ) *(R0+SUMA) 

W ( I ) = ( REAL ( IMPED ( I ) ) +SUMA+RO ) **2+ (AIMAG (IMPED ( I) )+SUMB) **2 

IF  (W(I)  .EQ.  0.0)  THEN 

T(I)=0.0 
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GO  TO  100 
ENDIF 

T(I)=TT(I)/W(I) 
100  CONTINUE 

DO  200  1=1, M 
F(I)=TPG-T(I) 
200  CONTINUE 
RETURN 
END 


C  *****************»*ratIONAL  FUNCTION**************** 

THIS  PROGRAM  DETERMINES  THE  RATIONAL  FUNCTION  OF  REISTIVE 
VALUES . 

EXTERNAL  EXAMP1 

INTEGER  M , N , IXJAC , NS IG , MAXFN , IOPT , INFER , IER, 

+1, J,K,L,NK,KK 

REAL  PARM (4 ) ,X(9)  ,F(11) ,XJAC(11,9) ,  XJTJ(45) ,  WORK(120) ,  EPS, 
+ DELTA,  SSQ, W (20) ,B(20) ,RQ(20) ,RX(20) ,A(1) 

COMMON  RQ,RX,W 

OPEN (UNIT=1 , FILE= ' RQM  DAT ' , STATUS=  * OLD ' ) 

OPEN (UNIT=2 , FILE* 1 FREQM  DAT • , STATUS® • OLD ' ) 

OPEN (UNIT=3 , FILE* ' RAT  DAT  B ', STATUS* *  OLD • ) 

M=ll 
N=9 

IXJAC=M 
NSIG=3 
EPS=0 . 0 
MAXFN=1000 
IOPT=l 
DELTA* 0 . 0 

DO  10  1*1, N 
10  X (I) *1. 0 

READ(1,*) (RQ ( I ) , 1=1 , M) 

READ ( 2 , * ) (W(I) ,1=1, M) 

CALL  XSSQ ( EXAMP1 , M , N , NS IG , EPS , DELTA , MAXFN , IOPT , PARM , X , SSQ , F , 
rXJAC, IXJAC, XJTJ, WORK, INFER, IER) 

L=(N-l)/2 
A ( 1) =X (1) 

B(1)=X(2) **2+2* (X(3)) 

WRITE (3,*) 'A(l) ' ,A(1) , 'B(l) ' ,B(1) 

DO  20  1*3, L 
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SUM=0 . 0 
K-I 

PRINT  * ,  K 
DO  30  J=3,K 

SUM=SUM+X ( J-l) *X(2*K-J+1) 

30  CONTINUE 

B(I)=X(I)**2  +  (2* (X(2*K-1)  +  SUM)) 

WRITE (3 ,  *) 'B(I) ' , B (I) 

20  CONTINUE 

B(4)=X(N)**2 

WRITE ( 3  ,  * ) '  B(4) » ,B(4) 

WRITE (3,*) 'X'  ,X 
WRITE (3 , *) 1 SSQ  * , SSQ 
WRITE (3,*) 'RQ* ,RQ 
WRITE ( 3  ,  * ) ' RX ' , RX 
END 

C  ***********  CALLING  SUBROUTINE  FOR  RATIONAL  FUNCTION  ******** 
SUBROUTINE  EXAMP1 (X, M, N, F) 

INTEGER  M, N, I , J 

REAL  X (N) , F (M) , SUMA, SUMB, SUMC, RX ( 20 ) ,RQ(20) ,W(20) 

COMMON  RQ , RX , W 

DO  10  1=1, M 
SUMA=0 . 0 
SUMB=0 . 0 
SUMC=0. 0 

DO  5  J=1 , N-l 

5  SUMA=SUMA+(X(J)**2)*(W(I)**2) 

SUMA=(X(1)**2) *(W{I) **2 ) 

DO  6  J=2,N 

SUMB=SUMB+X ( J) * (W ( I ) ** ( J-l ) ) 

6  SUMC=SUMC+X(J)*((-W(I))**(J-1)) 

RX ( I) =SUMA/ ( . 5* ( ( ( 1+SUMB) **2 ) + ( ( 1+SUMC) **2 ) ) ) 

10  CONTINUE 


DO  40  1=1, M 
F ( I ) =ABS (RX ( I ) -RQ ( I ) ) 
40  CONTINUE 
RETURN 
END 
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APPENDIX  B 
iriSL  SUBROUTINE 


FILE:  ZXSSQ  FORTRAN  A 


IriSL  ROUTINE  NAME  -  ZXSSQ 


COMPUTER 
LATEST  REVISION 
PURPOSE 

USAGE 

ARGUMENTS  FUNC 


n 

N 

NS  IG 


EPS 


DELTA 


MAXFN 

10  PT 


PARM 


-  I BM/S INGLE 

-  NOVEMDER  1,  1904 

-  MINIMUM  OF  THE  SUM  OF  SQUAHES  OF  M  FUNCTIONS 

IN  N  VARIABLES  USING  A  FINITE  DIFFERENCE 
LEVENEEBG-MARQUARDT  ALGORITHM 

-  CALL  ZXSSQ(PUNC,H.N,NSIG, BPS. DELTA, MAXFN, IOPT, 

PARH.X, S3Q,  P,  XJAC,IXJAC,XJTJ,  WC KK  ,  1 II F  ER ,  I  ER  j 


A  USES  SUPPLIED  SUBROUTINE  WHICH  CALCULATES 
THE  RESIDUAL  VECTOR  P { 1 J  ,  F (2)  ,  .  .  .  ,  P  (M)  FOR 
GIVEN  PARAMETER  VALUES  I  (  1)  ,X  (2)  ,  .  .  .  ,  X  (N) 

THE  CALLING  SEQUENCE  UAS  T  - 

CALL  FUNC <X,M.N . PI 
WHERE  X  IS  A  VECTOR  OF 
A  VECTOR  OF  LENGTH  M 
FUNC  MUST  APPEAR  III  AN 
IN  1HE  CALLING  PROGRAM. 

FUNC  MUST  NOT  ALTER  THE  VALUES  OF 
X  II),  1=1, ..  .  ,  N ,  H,  OR  N. 

THE  NUMBER  OF  RESIDUALS  OP  OBSERVATIONS 
^ I N  PUT) 


{i§  ^oltowiRo  fOfm 

LENGTH  N  AND  F  IS 
EXTERNAL  STATEMENT 


(INPUT) 


JT) 

IP  OH 


TWO 


THE  NUMBER  OP  UNKNOWN  PARAMETERS  M 

FIRST  CONVERGENCE  CRITERION.  IINPU1 
CONVERGENCE  CONDITION  SATISFIED  II 
SUCCESSIVE  ITERATIONS,  THE  PARAMETER 
ESTIMATES  AGREE,  COMPONENT  BY  COMPONENT, 

TO  NSIG  DIGITS. 

SECOND  CONVERGENCE  CRITERION.  (INPUT) 

CONVERGENCE  CONDITION  SATISFIED  IF,  ON  TWO 
SUCCESSIVE  ITERATIONS  THE  RESIDUAL  SUM 
OF  SQUARES  ESTIMATES  HAVE  RELATIVE 
DIFFERENCE  LESS  THAN  OR  EQUAL  TC  EPS.  EPS 
HAY  BE  SET  TC  ZERO. 

Till P D  CONVERGENCE  CRITERION.  (INPUT) 

CONVERGENCE  CONDITION  SATISFIED  IP  THE 
(EUCLIDEAN)  NORM  CF  THE  APPROXIMATE 
GPADIENT  IS  LESS  THAN  OR  EQUAL  TO  DELTA. 
DELTA  HAY  DE  SET  TC  ZERO. 

NOTE,  THE  ITERATION  IS  TERMINATED.  AND 
CONVERGENCE  IS  CONSIDERED  ACHIEVED,  IF 
ANY  ONE  OP  THE  THREE  CONDITIONS  IS 
SATISFIED. 

INPUT  MAXIMUM  NUMBER  OP  FUNCTION  EVALUATIONS 
(I. E.,  CALLS  TO  SUBROUTINE  PUNC)  ALLOWED. 

THE  ACTUAL  NUMBER  OF  CALLS  10  FUNC  MAY 
EXCEED  HAXPN  SLIGHTLY. 

INPUT  OPTIONS  PARAMETER. 

IOPT  =  0  IMPLIES  BROWN'S  ALGORITHM  WITHOUT 
STRICT  DESCENT  IS  DESIRED. 

I0PT=1  IMPLIES  STRICT  DESCENT  AND  DEFAULT 
VALUES  FOR  INPUT  VECTOR  PARM  ARE  DESIRED. 

I0PT=2  IMPLIES  STRICT  DESCENT  IS  DESIRED  WITH 
USER  PARAMETER  CHOICES  IN  INPUT  VECTOR  PABH. 

INPUT  VECTOR  OF  LENGTH  4  USED  ONLY  FOR 

10 PT  EQUAL  TWO.  PARM  (I)  CONTAINS.  WHEN 
1=1.  THE  INITIAL  VALUE  OF  THE  NARQUAROT 
PARAMETER  USED  TO  SCALE  THE  DIAGONAL  OF 
THE  APPROXIMATE  HESSIAN  MATRIX,  XJT.T. 

UY  THE  FACTOR  (1.0  ♦  PARfl(l)).  A  SMALL 
VALUE  GIVES  A  NEWTON  STEP.  WHILE  A  LARGE 
VALUE  GIVES  A  STEEPEST  DESCENT  STEP. 

THE  DEFAULT  VALUE  POR  PARN(1)  IS  0.01. 

1*2,  THE  SCALING  FACTOR  USED  TO  NODIPY  THE 
MARQUARDT  PARAMETER,  WHICH  IS  DECREASED 
BY  PARM  (2)  AFTER  AN  IMMEDIATELY  SUCCESSPUI 


zxssggio 

ZX.5SGG20 
ZXSSOO JO 
ZXJSQ&40 
ZXSSOOs J 
ZXSSOGou 
ZXSS0070 
ZXSSCOUO 
ZXSS0090 
ZXSSO  100 

zxsso  no 

ZXSSGIZO 
2XSS01  JO 
ZXSS0140 
ZXSSO  150 
ZXSSO  160 
ZXSSO  170 
ZXSSO  100 
ZXSSO  190 
ZXSS0200 
ZXSS02  1C 
ZXSSO  220 
ZXSSO  2 JO 
ZXSS0240 
ZXSS025C 
ZXSSJ2oO 
ZXSSO 2 70 
ZXSS02H0 
ZXSS0290 
ZXSSO  30 C 
ZXSS0310 
ZXSSO  320 
Z  XS  SO  3  30 
ZXSSO  3«0 
ZXSSO  350 
ZXSSO  3o 0 
ZX3S0370 
ZXSSO  3<J0 
ZXSSO  390 
ZXSS04J0 
ZXSS0410 
ZXSS0420 
ZXSS04  30 
ZX3S0440 
ZXSS0450 
ZXSS04C0 
ZXSS0470 
ZXSS04d0 
2XSS0490 
ZXSS0500 
ZXSS0S 1 0 
ZXSS0520 
ZXSS0530 
ZXSS0540 
ZXSS0550 
ZXSS05o0 
ZXSSG570 
ZXSS05U0 
ZXSSO  5  90 
Z  XS  SO  bO  0 
ZXSS06  10 
2XSS0620 
Z  XS  SOb  30 
ZX3S0640 
ZXSS0650 
ZXSS0660 
Z  XS  SO  6  70 
ZXSS06B0 
ZXSS0690 
ZXSS0700 
ZXSS07  1C 
ZXSS0720 
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KILE:  ZXSSQ 


FORTRAN  A 


X 

ssq 

F 

XJ  AC 
IX  J  AC 

XJTJ 

WORK 


INFER 


IER 


DESCENT  DIRECTION,  AND  INCREASED  BY  THE 
SQUARE  OF  PARM  (2)  IF  NOT.  PARM  (2)  BUST 
HE  GREATER  THAN  ONE,  AND  TWO  IS  DEPAULI. 
1=3,  AN  UPPER  BOUND  FOR  INCREASING  THE 
fUKOUARDT  PAR  Art ETER .  THE  SEARCH  FOR  A 
DESCENT  POINT  IS  ABANDONED  IP  PARM  (3)  IS 
EXCEEDED.  PARMJ3)  GREATER  THAN  100.0  IS 
RECOMMENDED.  DEPAULT  IS  120.0. 

1=4,  VALUE  EOR  INDICATING  WHEN  CENTRAL 

RATHER  THAN  FORWARD  DIFFERENCING  IS  TO  BE 
USED  FOR  CALCULATING  THE  JACOBIAN.  THE 
SWITCH  IS  MADE  WHEN  THE  NORM  OF  THE 
GRADIENT  OF  TUB  SUM  OF  SQUARES  FUNCTION 
BECOMES  SMALLER  THAN  PARM  (4) .  CENTRAL 
DIFFERENCING  IS  GOOD  IN  THE  VICINITY 
OF  THE  SOLUTION,  SO  PARM  (4)  SHOULD  BF. 
SMALL.  THE  DEPAULT  VALUE  IS  0.10. 

VECTOR  OF  LENGTH  N  CONTAINING  PARAMETER 
VALUES. 

Oil  INPUT,  X  SHOULD  CONTAIN  THE  INITIAL 
ESTIMATE  OF  THE  LOCATION  OF  THE  MINIMUM. 

ON  OUTPUT,  X  CONTAINS  THE  PINAL  ESTIMATE 
OF  THE  LOCATION  OF  THE  MINIMUM. 

OUTPUT  SCALAR  WHICH  IS  SET  TO  THE  RESIDUAL 
SUMS  OF  SQUARES,  F (1)  ♦* 2+ . . . ♦ P  (M) *♦ 2 ,  FOR 
THE  FINAL  PARAMETER  ESTIMATES. 

OUTPUT  VECTOR  OF  LENGTH  M  CONTAINING  THE 

RESIDUALS  FOR  THE  FINAL  PARAMETER  ESTIMATES 

OUTPUT  M  BY  N  MATRIX  CONTAINING  THE 

APPROXIMATE  JACOBIAN  AT  THE  OUTPUT  VECTOR  X 

INPUT  ROW  DIMENSION  OF  MATRIX  XJAC  EXACTLY 
AS  SPECIFIED  IN  THE  DIMENSION  STATEMENT 
IN  THE  CALLING  PROGRAM. 

OUTPUT  VECTOR  OF  LENGTH  [N»1)*N/2  CONTAINING 
THE  N  BY  N  MATRIX  ( XJ AC-T fi A  NS  PC SE D)  *  (XJAC 
IN  SYMMETRIC  STORAGE  MODE. 

WORK  VECTOR  OF  LENGTH  5*N  ♦  2* M  ♦  (N»1)*N/2. 

ON  OUTPUT,  WORK  |I)  CONTAINS  FOR 

1  =  1,  THE  NO  BN  OF  THE  GRADIENT  DESCRIBED 
UNDER  INPUT  PARAMETERS  DELTA  AND  PAFM  (4) . 
1=2,  THE  NUMBER  OF  FUNCTION  EVALUATIONS 
REQUIRED  DURING  THE  WOhK(5)  ITERATIONS. 
1=3,  THE  ESTIMATED  NUMBER  OF  SIGNIFICANT 
DIGITS  IN  OUTPUT  VECTOR  X. 

1=4,  THE  FINAL  VALUE  OF  THE  MARQIJARDT 

SCALING  PARAMETER  DESCRIBED  UNDER  ?ARM(1) 
1=5,  THE  NUMBER  OF  ITERATIONS  (I.E.,  CHANGE 
TO  THE  X  VECTOR)  PERFORMED. 

SEE  PROGRAMMING  NOTES  FOR  DESCRIPTION  OF 
THE  LATTER  ELEMENTS  OF  WORK. 

AN  INTEGER  THAT  IS  SET,  ON  OUTPUT,  TO 

INDICATE  WHICH  CONVERGENCE  CRITERION  WAS 
SATISFIED. 

INFER  =  0  INDICATES  THAT  CONVERGENCE  FAILED. 
IER  GIVES  FURTHER  EXPLANATION. 

INFER  =  1  INDICATES  THAT  THE  FIRST  CRITERION 


WAS  SATISFIED 
INFER  =  2  INDICATES  THAT  THE  SECOND  CRITERION 
WAS  SATISFIED. 

INFER  =  4  INDICATES  THAT  THE  THIRD  CRITERION 
WAS  SATISFIED. 

IP  MORE  THAN  ONE  OP  THE  CONVERGENCE  CRITERIA 
WERE  SATISFIED  ON  THE  FINAL  ITERATION, 

INFER  CONTAINS  THE  CORRESPONDING  SUM. 

( E. G.  ,  INFER  =  3  IMPLIES  FIRST  AND  SECOND 
CRITERIA  SATISFIED  SIMULTANEOUSLY). 

ERROR  PARAMETER  (OUTPUT) 

TERMINAL  ERROR 

IER=129  IMPLIES  A  SINGULARITY  WAS  DETECTED 
IN  THE  JACOBIAN  AND  RECOVERY  FAILED. 

IBR= 1 30  IMPLIES  AT  LEAST  ONE  OF  N,  N,  IOPT, 
P ARM  ( 1)  ,  OR  PARM  (2)  WAS  SPECIFIED 


ZXSSQ730 
2XSS074Q 
2XSS0750 
ZXS  30  ?b  0 
ZXSS077  0 
ZXSS0780 
ZXSS0790 

zxssoaoo 

ZX3S0810 
ZXSS0820 
ZXSS0630 
ZXSS084 0 
ZXSS0850 
ZXSS0860 
2XSS0870 
ZXSS088Q 
ZXSS0890 
ZXSS0900 
ZXSS0910 
ZXS  SO  9  20 
ZXSS09  30 
ZXSS0940 
2X5S0950 
2XSS0960 
ZXSS0970 
ZXSS0980 
ZXSS0990 
. ZXSS  1000 
ZXSS  1010 
.ZXSS1020 
ZXSS1030 
ZXSS1040 
ZXSS  1 0  50 
ZXSSIOoO 
) ZXSS1070 
ZXSS10B0 
ZXSS1090 
ZX5S1 100 
ZXSS1  1  10 
ZXSS1  IzO 
ZX3S11  30 
ZXSS1140 
ZXSS1 150 
ZXSS1  loO 
ZXSS  1170 
.ZXSS1  100 
SZXSS1190 
ZXSS  1 200 
ZXSS1210 
ZXSS1220 
ZXSS 1230 
ZXSS  1  240 
ZXSS 1 250 
ZXSS 1 260 
ZXSS1270 
ZXSS 1280 
ZXSS 1290 
ZXSS1300 
ZXSS13  10 
ZXSS 1 320 
ZXSS 1 330 
2XSS1  340 
ZXSS  1350 
ZXSS136C 
ZXSS1370 
ZXSS  1  380 
ZXSS  1 390 
ZXSS  1400 
ZXSS  1  4  1C 
ZXSS1420 
ZXSS1430 
ZXSS  1440 
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PILE:  ZXSSQ 


POSTS*  li  A 


INCORRECTLY . 
I£R=132  IMPLIES 
RECOVERY  FROM 
VECTCf.  X  HAS 
FIRST  SINGULA 
IER=133  IMPLIES 
WARNING  ERROR 
IER  =  38  IMPLIES 
THE  SOLUTION 
I ER=  39  IMPLIES 
PARAMETER  EXC 
USUALLY  MEANS 
ACCURACY  WAS 


THAT  AFTER  A  SUCCESSFUL 
A  SINGULAR  JACOBIAN,  THE 
CYCLED  BACK  TO  THE 
BUY. 

THAT  SAXFN  WAS  EXCEEDED. 

THAT  THE  JACOBIAN  IS  ZERO. 
X  IS  A  STATIONARY  POINT. 
THAT  THE  HARQUARDT 
EEDED  PARMJ3).  THIS 
THAT  THE  REQUESTED 
NOT  ACHIEVED. 


PRECISION/HARDWARE  - 


SINGLE  AND  DODBLE/B32 
SINGLE/H36,H48,U60 


REQD.  I  MS  L  ROUTINES  -  L  EQT 1  P,  Lfl  DECP,  LUE  LMP,  UERSET  ,  UERTST  ,  UG  ET 10 


NOTATION 


COPYRIGHT 

WARRANTY 


INFORMATION  ON  SPECIAL  NOTATION  AND 

CONVENTIONS  IS  AVAILABLE  IN  THE  MANUAL 
INTRODUCTION  OB  THROUGH  IS 5L  ROUTINE  QHELP 

1982  BY  IMSL,  INC.  ALL  RIGHTS  RESERVED. 

IMSL  WARRANTS  ONLY  THAT  IMSL  TESTING  HAS  BEEN 
APPLIED  TO  THIS  CODE.  NC  CTHER  WARRANTY, 
EXPRESSED  08  IMPLIED,  IS  APPLICABLE. 


SUBROUTINE  ZXSSQ 


INTEGER 

REAL 


INTEGER 


(FUNC,H,N,NSIG, EPS, DELTA, M A XPN, IOPT , P AR M , 

X, SSQ, f,XJ AC, IXJAC.XJTJ, WORK, INFER, IER) 
SPECIFICATIONS  FOR  ARGUMENTS 
M,N,NSIG,MAXFN,I0P1, I X J A C, I N PER , IER 
EPS. DELTA. PARH  (1)  ,X(N)  ,SSQ,F  (M)  ,XJAC(1)  , 

XJT5  (1) ,w6rK  (1) 

X J  AC  USED  INTERNALLY  IN  PACKED  FORM 
SPECIFICATIONS  FOR  LOCAL  VARIABLES 
IHJC.IGRADI.IGRADL.IGRADU.IDELXI.IDELXL, 
IDELXU,ISCAL1, ISC ALL. ISCALU. I X  NEW  1 , I KN E WL , 
IXBAD1, IFPI1,IfPL.IFPU,IFMLl .IFML. IEVAL, 
XBAD,ISW,ITER. J.IO AC,I,K,L,I5, JS,LI,LJ, ICOUNT, 
IZERO, LEVEL, LEVOLD 
AL.C0NS2.DN0RH,DSQ, 

ERL2, ERL2X.F0, rOSQ . F0SQS4 , G, HALF, 
HH.ONE.ONEflO.CNEPS. CNESFO.AX, 

PREC.REL.  RHII,SIG,SQfHF,SSQCLD,SUM,  TEN, 

TENTH, XDIF.X HOLD, UP, ZERO, 
IDABS,BELCON,P01,TWu,HUNTW,DELTA2 
SIG/6. 3/ 

AX/O. 1/ 

PO 1, TENTH, H ALF. ZERO. CNE. ON EP5, TWO, 

TEN, HU  NT W, ON EP 10/0. 01,0. 1. 0.5,0. 0, 

1. ,1.5,2., 10.0.1.2E2.1.E10/ 

ERROR  CHECKS 

FIRST  EXECUTABLE  STATEMENT 


I  EH  =  0 
LEVEL  =  0 

CALL  UERSET  (LEVEL, LEVOLD) 

IF  (M.LE. O.OR.M.GT. IX JAC . OR . N . LS. 0. OR .ICPT.  IT.  0 .  OR . IOPT .  GT.  2) 

*  OO  TO  305 
IMJC  =  IXJAC-,1 
IF  (IOPT. NE . 2)  GO  TO  5 

IF  (PARM  (2)  .  I.  E.  ONE.  OB.PARM  (1)  .  L  E.  ZERO)  GO  TO  305 

MACHINE  DEPENDENT  CONSTANTS 

>  PREC  »  TEN**  (-S1G-GNE) 

PEL  =  TEN**  <-SIG*HALF) 

RELCON  =  TEN*  * ( -N  SIG) 

WORK  VECTOR  IS  CONCATENATION  OF 
SCALED  HESSIAN, GRADIENT, DEL X, SC ALE, 
XNEW , XBA  D, F  (X  +  DEL)  ,F(X-DEL) 


))  GO  TO  305 
DEPENDENT  CONSTANTS 


ZXSS1450 
ZXSS1460 
ZXSS1470 
ZXSS1480 
ZXSS1490 
2XSS1500 
ZXSS1510 
ZXSS152C 
ZXSS1530 
ZXSS 1540 
ZXSS1550 
ZXSS1560 
ZXSS1570 
ZXSS1580 
ZXSS 1 590 
ZXSS 1 600 
ZXSS  1  6  1  0 
ZXSS 1620 
ZXSS1630 
ZXSS164Q 
ZXSSIoSO 
ZXSSI60O 
ZXSS1670 
2XSS1680 
ZXSS1690 
ZXSS1700 
ZXSS1710 
ZXSS1720 
ZXSS17J0 
-ZXSS1740 
ZXSS 1750 
ZXSS  1760 
ZXSS1770 
ZXSS 1780 
ZXSS1790 
ZXSS1800 
ZXSS1810 
ZXSS1820 
ZXSS 1 830 
ZXSS 1 8  4  0 
ZXSS 1850 
ZXSS1860 
ZXSS1870 
ZXSS 1880 
ZXSS 1 890 
ZXS519Q0 
ZXSS  1 9  1 0 
ZXSS 1920 
ZXSS 1930 
ZXSS 1940 
ZXSS  1950 
ZXSS  1  9t>0 
ZXSS 1 970 
ZXSS 1 980 
ZXSS  1  990 
ZXSS2Q00 
ZXSS20 10 
ZXSS2020 
ZXSS20  3C 
ZXSS2040 
ZXSS20  50 
ZXSSZOoO 
ZXSS2070 
ZXSS208C 
ZXSS2090 
ZXSS2 100 
ZXSS2110 
ZXSS2 120 
ZXSS21 JO 
ZXSS2 1 40 
ZXSS2150 
ZXSS2160 


40 


FILE:  ZXSSQ 


POP TP  AN  A 


IGHAD1 
IGR ADL 
IGft  ADO 
IDELX  1 
IDELXL 
IDELXU 
ISCAL  1 
ISCALL 
ISCALU 
IXNEW  1 
1XNEWL 
IXBAD1 
IPPL1  = 


=  IUHAD « 

=  IG8ADUN 
=  IGft ADO 
=  IDELX 1 ♦ 1 
=  IDELX 1 »N 
=  IDELXU 
=  ISCALU  1 
=  ISCALUN 
=  ISCALU 
=  IXNEWU1 
=  IX  NEW  1 *N 
IXBADUN 


IPPL  =  IPPLU1 
IFPU  =  IFPLUrt 
IPML1  =  IPPU 
IPNL  =  IFNLU1 
IKJC  =  IXJAC  - 


INITIALIZE  VARIABLES 


IP 

AL 

PO 

UP 


GO 

GO 


TO 

TO 


2D 

10 


AL  =  ONE 
C0NS2  =  TENTH 
IF  (IOPT.EQ.0) 
llOPT.EQ. ij 

=  PARS  1) 

=  PARS  2 
=  PARK  3 
C0NS2  =  PARH(4) 

GO  TO  15 
10  AL  =  P01 
PO  =  TWO 
UP  =  HUNTW 
15  ONESPO  =  ONE/P0 
FOSQ  =  FO*PO 
F0SQ54  *  P0SQ**4 
20  IEVAL  *  0 

DELTA  2  =  DELTA*HA  LF 

EBL2  =  ONEP10 

IBAD  =  -99 

ISW  =  1 

ITER  =  -1 

INFER  =  0 

IER  =  0 

DO  25  J=IUELXL, IDELXU 
WORK  (J)  =  ZERO 

25  CONTINUE 
GO  TO  165 

MAIN  LOOP 

30  SSQOLD  =  S3Q 

CALCULATE  JACOBIAN 

IP  (INPER.GT.O.OR. IJAC.GE.N.OH. 10 PT . Eg. 0 . 0 R . ICOUNT . GT . 0) 

RANK  ONE  UPDATE  TO  JACOBIAN 

I JAC  =  IJACH 
DSQ  =  ZERO 

DO  35  J=IDELXL. IDELXU 

DSQ  =  DSQ*WOBK  (J)  *WORK  (J) 

35  CONTINUE 

IF  (DSQ . L E . ZERO)  GO  TO  55 
DO  50  1=1, H 

G  =  F(l)  -WORK  (IFNT.UI) 

K  =  I 

DO  40  J=I DEL  XL. IDELXU 
G  =  G+XJAC (K) AWORK (JJ 
K  =  Kt IXJAC 
40  CONTINUE 

G  =  G/DSQ 
K  =  I 

DO  45  J=IDELXL, IDELXU 

X  JAC  (K1  =  XJAC(K)  -G*UORK  (J) 

K  =  K  *1 X  J  AC 
45  CONTINUE 

50  CONTINUE 
GO  TO  ao 


GO  TO  55 


ZXSS2170 
ZXSS2100 
ZX3S2190 
ZXSS227G 
ZXSS2210 
ZXSS22^:0 
ZXSS2230 
ZXSS2240 
ZXSS2250 
Z  XSS  2  2o  0 
ZXSS2270 
ZXSS22«0 
ZXSS2290 
ZXSS2300 
ZXSS2310 
ZXSS232C 
ZXSS2330 
ZXSS2340 
ZXSS2350 
ZXSS2360 
ZXSS2370 
ZXSS238G 
ZXSS2390 
ZXSS2400 
ZXSS24  10 
ZXSS24  20 
ZXSS2430 
ZXSS2440 
ZXSS2450 
ZXSS24bC 
ZXSS24  70 
Z  XSS  248  0 
ZXSS2490 
ZXSS2500 
ZXSS2510 
Z  XSS2520 
ZXSS2530 
ZXSS254  0 
ZXSS255C 
ZXSS2560 
ZXSS2570 
ZXSS2580 
ZXSS259C 
ZXSS2600 
ZXSS2610 
ZXS52620 
ZXSS2630 
ZXSS2640 
ZXSS2650 
ZXSS2bo0 
ZXSS2670 
ZXSS2680 
ZXSS2690 
ZXSS2700 
ZXSS2710 
ZXSS2720 
ZXSS2730 
ZXSS2740 
ZXSS2750 
Z  XSS276  0 
ZXSS2770 
ZX5S27B0 
ZXSS2790 
ZXSS2U00 
ZXSS2310 
ZXSS2020 
ZXSS2830 
ZXSS2840 
Z  XS  S  2850 
ZXSS2860 
Z XSS 28 70 
ZXSS2880 


41 


i  IL£:  ZXSSQ 


FORTRAN  A 


5-)  I JAC  =  0 
K  =  -  III  JO 
DO  75  J=1.N 
K  =  KttM.JC 
XDAUS  =  AH  5 { X 1 J)  1 
HH  =  REL*{ArtAX1 (XDAUS, AX) ) 
XiiOLD  -  X{J) 

X  (J)  =  X  (J)  *  HH 

CALL  FUNC  (X.M.N .WORK  (IFPL) ) 
IEVAL  =  IEVAtn 


J A CO  El  AN  BX  INCREMENTING  X 


X(J)  =  XHOLD 
IF  (ISU. EQ. I)  GO  TO  65 

X(J)  =  XHOLD-HH 


CENTRAL  DIFFERENCES 


60 

65 

70 


ALL  FUNC  {X,M.N,WOIiK(IFHL}  ) 

IEVAL  =  IEVAL+1 
X{J)  =  XIIOLD 
RHH  =  HALF/IIH 
DO  60  L=IFPL,IFPU 
K  =  K  ♦  1 

XJACfK)  =  (WORK  (I) -WORK  (I»M)  J  *HIIH 
CONTINUE 
GO  TO  75 


75  CONTINUE 
HO 


RHH  =  ONE/Htl 
DO  70  1=1. M 
K  =  K  ♦  l 
XJACfK)  = 
CONTINUE 


ERL2X  =  ESL2 
ERL  2  =  ZERO 
K  =  -IMJC 

DO  90  J=IGRADL,  IGRAOU 
K  =  K » IMJC 
SUM  =  ZERO 
DO  85  1=  1 .  M 
K  =  K  ♦  i 

SUM  =  SUH+XJAC  (K)  *F  (I) 
85  CONTINUE 

WORK  (J)  =  SUM 

ERL2  =  ERL2«-SUM*SUM 
90  CONTINUE 

ERL2  =  SORT  (ERL2) 


FORWARD  DIFFERENCES 

(WORK  (LFPLWI)-P  (I)  )  * Rllll 

CALCULATE  GRADIENT 


IF  (IJAC.GT.O)  GO  TO  95 
IF  (ERL  2.  I.E  .  DELTA  2)  INFER  =  1NPER*4 
IP  (ERL2. LE.CCNS2) 


CONVERGENCE  TEST  POL  NORM  CF 


ISW  =  2 


95  L  = 
IS 
DO 


CALCULATE  THE  LOWER  SUPER  TP 
JACOBIAN  (TRANSPOSED)  *  JACO 


JS 

DO 


CONTItU 


100 
105 

1  17  CONTINUE 

IP 
IF 


=  -IX  J  AC- 
110  1=1,  N 
IS  =  I3HXJAC 
=  - IX  J  AC 
105  J=1,I 
JS  =  1S*I XJ  AC 
L  =  L*1 
SUM  =  ZERO 
DO  100  K=1,N 
U  =  ISM 
LJ  =  JS»K 

SUM  =  SUM»XJAC (LI) *XJAC(LJ) 
CONTINUE 
X.ITJ^L)  =  SUM 


IF  (IN 

(IEVAL. GE.  MAXFN) 


FER.GT.O)  GO  TO  315 
-■*  .  GO  TO  290 


CONVERGENCE  CHECKS 


ZXSS2890 
ZXSS2900 
2XSS29 10 
ZXSS2920 
ZXSS2930 
ZXSS29N0 
ZXS32950 
ZXSS2960 
ZXSS2970 
2X5S29U0 
ZXSS2990 
ZXSS3000 
ZXSS30 10 
ZXSS3020 
ZXSS30 JO 
2XSS304O 
ZXSS3050 
ZXSS3060 
ZXSS3070 
ZXSS3080 
ZXSS3090 
2XSS3100 
ZXSS3  11C 
ZXSS3 120 
ZXSS313C 
ZXSS3140 
ZXSS3 150 
ZXSS3 160 
ZXS53170 
ZX3S3180 
ZXSS3190 
ZXS5320C 
ZXSS32 1 0 
ZXSS3220 
ZXSS3230 
ZXSSJ24C 
ZXSS3250 
ZXSS32oO 
ZXS53270 
ZXSSJ280 
ZXSS3290 
ZXSS3300 
ZXSS3310 
ZXSS3320 
ZXSS3330 
ZXSS3340 
GRADIENT ZXSSJ35Q 
ZXSSJ360 
ZXSS3370 
ZXSSJ3H0 
IANUE  OF  ZXSS3390 
81  AH  ZXSS3400 

ZXSSJ4  10 
ZXSS3420 
ZXSS34 30 
ZXSS3440 
ZXSS 3450 
ZXSS34uO 
ZXSS3470 
ZXSS J4d0 
ZXSS3490 
2XSS3500 
ZXSS3510 
ZXSS 3520 
ZXSS3530 
ZXSS3540 
ZXSS3550 
ZXSSJSoO 
ZXSS3570 
Z  XS  S  3580 
ZXSS3590 
ZXSS3600 


FILE;  ZXSSQ 


FORTRAN  A 


COMPOTE  SCALING  VECTOR 

IF  ( IOPT . EQ . 0)  GO  TO  120 

DO  115  .1=1,  N 
K  =  K*J 

WORK  (ISCAL  1  t.J)  =  XJT.I(K) 

CONTINUE 
GO  TO  135 

CONFUTE  SCALING  VECTOR  AND  NOUN 

DNO  BN  =  ZEhO 
K  =  0 

DO  125  J= 1 , N 
K  =  K+J 

WORK  {ISC  'L  1  ♦  J)  =  SORT  ( X  JT.I  (K)  ) 

DNORi  =  DNCRM*XJTJ  £K)  wXJTJ  K) 

CONTINUE 

DNO  Bfl  =  ONE/SQRT  (DNORN) 

NORMALIZE  SCALING  VECTOR 

DO  130  J=ISCALL, ISCAL U 

WORK(J)  =  WORK  (J)  *  DNOiifl*  ER  L2 
CONTINUE 

ADD  L-M  FACTOR  TO  DIAGONAL 

ICOUNT  =  0 
K  =  0 

DO  150  1=1, N 
DO  145  .1  =  1,1 
K  =  K  *  1 

WORK  { K)  =  XJTJ  (K) 

CONTINUE 

WORK(K)  =  WORK  (K)  »WORK  (ISCALUI)  *AI 
WORK  [I DELX 1 ♦ I)  =  WORK { IG  RADI ♦  I ) 

CONTINUE 

CiiCLESKY  DECOMPOSITION 
CALL  LEQT1P  ( WORK , 1 , N . WOR K ( IDEL XL)  , N , 0 , G , XHOLD, I E R) 


CALCULATE  SUN  OF  SQUARES 


IF 

(IER . EQ. 0) 

GC 

TO 

160 

IER 

=  0 

IF 

(I J AC. GT.0) 

GO 

TO 

55 

IP 

(IBAD.  LE  .0) 

GO 

TO 

240 

IP 

(IBA0.GE.2) 

GO 

TO 

310 

GO 

TO  190 

IF  (IBAD.NE.-9U)  IBAD  =  0 

CALCULATE 

DO  170  J  = 1 , H 

WORK  (IXNEWUJ)  =  X  (J) -WORK  (IDELXUJ) 
CONTINUE 

CALL  FUNC  (WORK  (IXNBWL)  ,  (1 ,  N  ,  WCfi  K  ( I F  PL) ) 

IEVAL  =  I  £ V  AL ♦  1 

SSQ  =  ZERO 

OO  175  1  =  IF  PL , IPPU 

SSQ  =  SSQ*- WORK  (I)  *WORK  (I) 

CONTINUE 

IF  (ITEH.GE.0)  GO  TO  185 


ITER  =  0 
SSQCI.D  =  SSQ 
DO  180  1=1, H 

F  ( 1)  =  WORK (IF  PL  1*1) 

CONTINUE 
GO  TO  55 

IF  (IOPT.EQ.0)  GO  TO  215 
IF  (SSQ.LE.  SSQCLD)  GO  TO  205 


SSQ  FOR  INITIAL  ESTIMATES  OF  X 


ICOUNT  =  ICOUNT* 1 
AL  =  A  I.  *  F  0  S  Q 

IF  t J  AC. EQ . 0)  GO  TO  195 

IP  ICOUNT. GE. 4. OR. AL.GT.UP)  GOTC  200 

IF  AL.LE.UP)  GC  TO  140 

IF  ISAD. EQ. 1)  GO  TO  310 

IEH  =  39 

GO  TO  315 

AL  =  AL/P0SQS4 


CHECK  DESCENT  PROPERTY 
INCREASE  "ARAMETER  AND  TRY  AGAIN 


ZXS  S  3t)  1  0 
ZX5S362C 
ZXSSJ630 
ZXSS3640 
ZXSS3650 
ZXSS  3 6b  C 
ZXSS  3b7  0 
ZXSS36B0 
ZXSS3690 
ZXSS  370  0 
ZXSSJ7 10 
ZXSS3720 
ZXSS3730 
ZXSS  37  4  0 
ZXSS375Q 
ZXSS3760 
ZXSS  37  70 
ZXSS370Q 
ZXSSJ7J0 
ZXSS3800 
ZXSS381 u 
ZXSSJ820 
ZXS3J830 
ZXSS3840 
ZXSS  3850 
zxssiaoo 
ZXSS3870 
ZXSSJ880 
ZXSS3890 
ZXSS3900 
ZXSS3910 
ZXSS3920 
ZXSS3930 
ZXSS3940 
ZXSS3950 
ZXSS3960 
ZXSS3970 
ZXSS3980 
ZXSS 3990 
ZXSS4000 
ZXSS4010 
ZXSS402C 
ZXS  S4  0  30 
ZXSS40  4  0 
ZXSS4050 
ZXSS406Q 
ZXSS4070 
ZXSS4080 
ZXSS4090 
ZXSS4 100 
ZXSS4  1  10 
ZXSS4 120 
ZXSS4  130 
ZXSS4 140 
ZXSS4  150 
ZXSSU  160 
ZXSS4 170 
ZXSS4 180 
ZXSS4 190 
ZXSS4200 
ZXSS42 10 
ZXSS4  220 
ZXSS42 JO 
ZXSS4240 
ZXSS4250 
ZXSS4260 
ZXSS427  0 
Z  XS  34  28  0 
ZXSS4290 
ZXSS4300 
ZXSS4  J  1  0 
ZXSS4320 


43 


FILE:  /iSSQ 


pohtban  A 


GO  TO  55 

ADJUST  (1 AHQOABDT  PAR AHETER 

IF  (ICOUNT.  SQ.  0)  AL  =  AL/FO 
IF  (ERL2X.LE. ZERO)  GO  TO  210 
G  *  ERL2/EEL2X 

IF  (EflL2.LT. ERL2X)  AL  =  AL* AH  AX  1  (CNESFO ,G) 

IP  (EEL  2. GT . EBL2X)  AL  =  A L* AH  IN  1  ( FO , G ) 

AL  =  AH AX  1  (AL, PBEC) 

ONE  ITERATION  CYCLE  COflPLEIED 

ITER  =  ITER ♦ 1 
DO  220  J= 1 , N 

X(J)  =  HORK(IXNEW1*J) 

CONTINUE 
DO  225  1*1, B 

WORK  (IFHL1+I)  =  F  ( I) 

Ffli  =  WORK  (IFPL 1*1 
CONTINUE 

RELATIVE  CONVERGENCE  TEST  FO?  X 

IF  JAI..GT.5.0)  GO  TO  30 
DO  2 JO  J=1,N 

XDIP  =  ABS(l/ORK(IDELX1+J))/Af1AX1  (ASS  (X  ( J )  )  ,  A  X) 

IF  (XDIF.G'f.  RELCOM)  GO  TC  235 
CONTINUE 
INFER  =  1 

RELATIVE  CONVERGENCE  TESI  FOR  SSQ 
SOD  IF  =  ASS  (SSQ-S'JQ0LD)/AHAX1  ISSQOl.D,  AX) 

I P  (SQDIF.LE.EPS)  INFER  *  INFEfi*2 
GO  1C  JO 

SINGULAR  DECOMPOSITION 

IF  (I QAD)  255,245 , 265 

CHECK  TO  SEE  IF  CURRENT 
ITERATE  RAS  CYCLED  BACK  TO 
THE  LAST  SINGULAR  POINT 

DO  250  J=1,N 

XIIOLD  =  WORK  (IXBADItJJ 

IF  (ABS  (X  (J) -XHOLD)  .GT. RELCON* AHAX 1  (AX, A BS  (XHOLD) ) )  GO  "O  255 
CONTINUE 
GO  TO  295 

UPDATE  THE  UAD  X  VALUES 

DO  260  J=1,N 

WORK  (IXBADUJ)  =  X  (J) 

CONTINUE 
IBAD  =  1 

INCREASE  DIAGONAL  CF  HESSIAN 

IF  (IOPT.NE.0)  GO  TO  2(30 
K  =  0 

DO  275  1*1, N 
DO  270  3*1, I 
K  *  K  ♦  1 

WORK  (  K)  *  XJTJ(K) 

CONTINUE 

WORK  (K)  =  ONEP5*  (XJTJ (K)  ♦  A L* ER 1 2*V CR K  (ISCAL1 ♦ I ) )  *REL 

CONTINUE 
IBAD  =  2 
GO  TO  155 

REPLACE  ZEROES  Otl  HESSIAN  DIAGONAL 

IZERO  =  3 

DO  205  J=ISCALL, ISCALU 

IF  (WORK  (J).GT.  ZERO)  GO  TO  285 
IZERO  *  IZEftCO 
WORK  (.1)  =  ONE 

CONTINUE 

IF  (IZERO. I.T.N)  GO  TO  143 
IER  =  38 
GO  10  315 

T ERN INAL  ER3CK 

IER  =  IEHH 

IER  ■=  I  E?  ♦  1 

IER  =  I  ER  ♦  1 

IER  =  I  ER  ♦  1 

IER  =  I ER ♦ 1 29 

IF  (IER.  EQ.  133)  GO  TO  335 


ZXSS4330 
ZXSS4 340 
ZXSS4350 
ZXSS4 360 
ZXSS437  3 
ZXSS4300 
ZXSS4390 
ZXSS4400 
ZXSS44 10 
ZXSS4420 
ZXSS44J0 
ZXSS4440 
ZXSS4450 
ZXSS4460 
ZXSS4470 
ZXSS44H0 
ZXSS4490 
ZXSS450G 
ZXSS45 10 
ZXSS4520 
ZXSS4530 
ZXSS4540 
ZXSS4550 
ZXSS4  5b0 
ZXSS4570 
ZXSS45B0 
ZXSS4593 
ZXSS4600 
ZXSS46  10 
ZXSS4b20 
ZXSS46 30 
ZXSS4640 
ZXSS4650 
ZXSS4660 
ZXSS4670 
ZXSS4b80 
2XSS4b90 
ZXSS47C0 
ZXSS47 10 
ZXSS4720 
ZXSS4730 
ZXSS4740 
ZXSS4750 
ZXSS4760 
ZXSS4770 
ZXSS4780 
ZXSS4790 
ZXSS4800 
ZXSS48  10 
ZXSS4820 
ZXSS48  30 
ZXSS4840 
ZXSS4850 
ZXSS4Ub0 
ZXSS4870 
ZXSS4880 
ZXSS4890 
ZXSS4900 
ZXSS4910 
ZXSS4920 
ZXSS49  30 
ZXSS4940 
ZXSS4950 
Z  XS  S4  9b  0 
ZXSS4970 
ZXSS4980 
ZXSS4  990 
ZXSS5000 
ZXSS5010 
ZX5S5020 
ZXSS5030 
ZXSS5040 


FILE:  ZXSSQ 


FORTRAN 


OUTPUT  ERL2,IEVAL,NSIC,AL, 

315  0  =  SIG 

00  I»OLD"''iBS(*OBK  l 

5F-lSSIi?^!;“:5fSSl,o,|ii.5fD..5Lo0,o(»»»xi|Sx,»BM«u1) 

J  20  CONTINUE 

I P  ( H  .  GT  .  2)  GO  TO  3  30 
DO  325  J  =  1,N 

325  WORK(J*5)  =  WORK  (JfIGRADI) 

330  WORK  (  1)  =  P.BL2*EuL2 
WORK  (2)  =  IEVAL 
SSO  3  SSQOLD 
WORK  (3)  =  G 
WORK  (4)  =  At 

WORK  (5)  *  ITER 

335  CALL  UERSKT  (LEVOLD,  LEVOI.D) 

IF  (IEE.EQ.  0)  GO  TO  9005 
9000  CONTINUE  , 

CALL  UERTST  (IgR,6HZX53Q  ) 

9005  RETURN 
END 


AND  ITER 


ZX3S5050 
ZXS35Q60 
ZX3S5070 
2  XSS50U  0 
2X335090 
ZX5S5100 
ZXSS5 110 
ZXSS5U0 
ZXSS5130 
ZXSS5140 
ZXSS5 150 
ZXSS5 16  0 
ZXSS5170 
ZXSS5  1U0 
ZXS35190 
ZX3S52QG 
ZX3352 1 0 
ZXSS522C 
ZXSS5230 
ZXSS5240 
ZX335250 
2  XS  S5  2b  0 


45 


LIST  OF  REFERENCES 


1.  Carlin,  Herbert  J.,  "A  New  Approach  to  Gain-Bandwidth 
Problems,"  IEEE  Transactions  on  Circuits  and  Systems,  Vol. 
CAS-24 ,  No.  4,  April  1977,  pp .  170-175. 

2.  Ramahi,  Omar  M.  and  Mittra,  Raj,  "Design  of  a  Matching  Network 
for  an  HF  Antenna  Using  the  Real  Frequency  Method,"  IEEE 
Transactions  on  Antennas  and  Propagation ,  Vol.  37,  No.  4, 
April  1989,  pp.  506-509 

3.  Proakis,  John  G. ,  Digital  Communications,  McGraw-Hill  Inc., 
1989,  pp.  151-152. 

4.  Bode,  H.W.,  Network  Analysis  and  Feedback  Amplifier  Design, 
New  York:  Van  Nostrand,  1945,  pp.  319. 

5.  Carlin,  Herbert  J.  and  Yarman,  Binboga  S.,  "The  Double 

Matching  Problem:  Analytic  and  Real  Frequency  Solutions," 

IEEE  Transactions  on  Circuits  and  Systems,  Vol.  CAS-3  0,  No.  1, 
January  1983,  pp.  22. 

6.  Kami,  Shlomo,  Network  Theory  and  Analysis  and  Synthesis , 
Allyn  and  Bacon  Inc.,  Boston,  1966,  pp.  190-193. 

7.  Beccari,  Claudio,  "Broadband  Matching  Using  the  Real  Frequency 
Technique,"  ISCAS  1984,  pp.  1231-1232. 

8.  Weeks,  Walter  L. ,  Antenna  Engineering,  McGraw-Hill  Book 
Company,  1968,  pp.  349-350. 

9.  Stutzman,  Warren  L.  and  Thiele,  Gary  A.,  Antenna  Theory  and 
Design,  John  Wiley  &  Sons,  1981,  pp.  310-313. 

10.  Rao,  B.  Rama,  "Optimized  Tapered  Resistivity  Profile  for 
Wideband  HF  Monopole  Antenna,"  IEEE  International  Antennas  and 
Propagation  Symposium,  1991  Digest,  Vol  2,  pp.  711-714. 

11.  Wu,  T.T.  and  King,  R.W.,  "The  Cylindrical  Antenna  with 
Nonreflecting  Restive  Loading,"  IEEE  Transactions  on  Antennas 
and  Propagation,  Vol.  AP-13,  May  1965,  pp.  369-373. 

12.  Stewart,  John  L. ,  Circuit  Theory  and  Design,  John  Wiley  and 
Sons,  Inc.,  New  York,  1956,  pp.  82-88. 


46 


13.  Janaswamy,  R.  Unpublished  Notes,  Naval  Postgraduate  School, 
Monterey,  California,  1992. 

14.  Davis,  W.  A.,  WIRE,  Virginia  Polytechnic  Institute  and  State 
University,  Blacksburg,  Virginia. 


47 


INITIAL  DISTRIBUTION  LIST 


Defense  Technical  Information  Center  2 

Cameron  Station 
Alexandria,  VA  22304-6145 

Library,  Code  52  2 

Naval  Postgraduate  School 
Monterey,  California  93943-5002 

Chairman,  Code  EC  1 

Department  of  Electrical  and 

Computer  Engineering 

Naval  Postgraduate  School 

Monterey,  California  93943-5000 

Prof.  R.  Janaswamy,  Code  EC/Js  1 

Department  of  Electrical  and 

Computer  Engineering 

Naval  Postgraduate  School 

Monterey,  California  93943-5000 

Prof.  David  C.  Jenn,  Code  EC/Jn  1 

Department  of  Electrical  and 

Computer  Engineering 

Naval  Postgraduate  School 

Monterey,  California  93943-5000 

COMNAVFORKOREA  1 

Unit  15250 

APO  AP  96205-0023 

Dr.  Felix  Schwering  1 

AMSEL-RD-C3 -  JA1 

Center  for  C3  Systems 

Fort  Monmouth,  New  Jersey  07703 


48 


